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ABSTRACT 

We have used the Green Bank Telescope to carry out a deep search for redshifted 
CO J=2-l line emission from an extended (>17 kpc) Lya blob (LAB), "Himiko", 
at z ~ 6.595. Our non-detection of CO J=2-l emission places the strong 3cr 
upper limit of L'^o < 1-8 x 10^° x (AF/250)i/2 ^ km s^^ pc^ on the CO line 
luminosity. This is comparable to the best current limits on the CO line luminos- 
ity in LABs at 2; ~ 3 and lower- luminosity Lyman-a emitters (LAEs) at z > 6.5. 
High-z LABs appear to have lower CO line luminosities than the host galaxies 
of luminous quasars and sub-mm galaxies at similar redshifts, despite their high 
stellar mass. Although the CO-to-H2 conversion factor is uncertain for galaxies 
in the early Universe, we assume Xco = 0.8 Mq (K km s~^ pc^)^^ to obtain the 
limit M(H2) < 1.4 x 10^*^ Mq on Himiko's molecular gas mass; this is a factor of 
> 2.5 lower than the stellar mass in the z ^ 6.595 LAB. 

Subject headings: galaxies: high-redshift, galaxies: ISM 



Introduction 



Significant populations of star-forming galaxies have recen tly been di scovered at high 



redshifts, z > 6, towards the end of the reionization epoch (e.g. lEUisI l2008l ). Many of these 



galaxies have been identified due to their excess eniission in narrow-band iniages centred on 



the r edshifted Lya wavelength (e.g. iRhoads et al.ll2000t iTaniguchi &: et al.ll2005l : IHu et al. 



2OIOI ): such systems are referred to as Lya emitters (hereafter LAEs). The star-formation 
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rates determined from measurements of the rest-frame UV continuum emission in LAEs at 

'1 



z ~ 6.5 tend to be low, SFR ~ 5 - 10 Mq yr 
of the emission region is typically small, ~ 1 kpc (ICowie et al.ll201ll ) 



Taniguchi &: et al.ll2005l ). and the size 



High-z LAEs thus appear to be relatively small galaxies, probably undergoing quiescent 
star formation. However, a recent Subaru narrow-band imaging survey has discovered a 
highly extended LAE at z ~ 6.6, with a spatial extent of > 10 kpc, far larger than that of 
the bulk of the z ~ 6.5 LAE population ( lOuchi et al.ll2009[ ). The object, "Himiko" , also has a 



far higher Lya line luminosity, > 10 erg s~ , than any other LAE at 2; > 6 . Both the large 
size and the high Lyman-a luminosity make Himiko similar to the "Lyman-g blobs" (LABs) 
that have so far been detected at z 2.5 — 4 (e.g. 



Dev et al. 


2005; 


Nilsson et al. 


2006: 



Steidel et al.ll2000l : iMatsuda et al.ll2004J : 
201l|). The LABs discovered at ^ ~ 3 



are extremely large (> 50 kpc in size) Lya-emitting nebulae that show a striking correlation 
with over-dense regions in the Universe. These have Lya lumin osities > 10"^^ erg s~ ^, similar 



to the luminosities seen in massive high-z radio galaxies (e.g. iReuland et al.ll2003l . but do 
not show radio emission. Follow-up multi-wavelength studies have shown tha t LABs tend 



Geach et al. 



Geach et al. 



to be associated with bright sub-millimeter or infra red-luminous galaxies (e.g. 
^05!) or obscured active galactic nuclei (AGNs; e.g. iBasu-Zych fc Scharj|2004j : 
2009,) . often with high star formation rates (~ 1000 M© per year). Recent Spitzer stndies 
have found that the infrared (IR) images of 60% of z ~ 3 LABs are consistent with an origin 
due to star formation, with the remaini ng 40% of the samp le likely to arise either due to 
AGN activity or an extreme starburst (IColbert et al.ll201l[). A few LABs have also been 



detected without a known associated bright galaxy (e.g. iNilsson et al.ll2006l ). 



No LABs have so far been detected at low redshifts, z <1 (IKeel et al.ll2009l ): the blobs 



thus appear to be a high-redshift phenomenon. The fact that LABs are typically located 



Steidel et al. 


2000; 


Erb et al. 


2011) 



(e.g. 

suggested for the sources that power the Lya emission (e.g. cooling inf lows of gas, outflows 
of ga s from starburst galaxies or AGNs, or eyen AGN photoionizat i on: iTaniguchi &: Shiova 



2000 



Haiman fc Reed I2OOII : iMatsuda et al.l 120041 : IPev et all 120051 : IPiikstra fc Loeb 



2009 



Geach et al.l l2009 l) , it appears unclear whether any single mechanism is capable of explain- 
ing the entire LAB population. However, while cooling radiation from gas strea ms may 
contribute some fraction of the LAB power (IHaiman et al.ll2000l : iGoerdt et al.ll2010l ). radia- 
tive transfer calcula tions have shown that such streams of halo gas cannot account for all 
the Lya luminosity ( iFaucher-Giguere et al.ll2010l ). Galaxy formation must hence be respon- 
sible for a significant fraction of the Lya emission from LABs. Studies of the star formation 
activity of galaxies in the vicinity of LABs are thus essential to understanding their nature. 
Observations of molecular and atomic gas tracers in the LABs would provide crucial insight 
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into the fuel for star formation. The most effective means of studying the cold molecu- 
lar g as reservoir in high- 2: galaxies is through observations of redshifted CO line emission 



[e.g. iSolomon fc Vanden Boutll2005l ). which provide a means of studying gas kinematics and 
estimating the total molecular gas mass available for star formation. 

At present, only two LABs have been searched for CO line emission, both at 2; ~ 3, but 
with no detections of molecular gas ( lYang et al.l 120121 ). We present here a search for CO 
J=2-l line emission in the z ~ 6.595 LAB, Himiko, in order to measure its cold molecular 
gas mass and constrain the kinematics of the interstellar gas in a massive galaxy during the 
epoch of reionizationl^ 



2. Himiko: a Lya blob at z=6.595 



The target of our CO line observations, "Himiko", is a giant 17 kpc) LAE at z = 
6.595, discovered in the SvLhain/ XMM-Newton Deep Survey ( lOuchi et al.ll2009l ). Himiko has 



an extremely high Lya line luminosity, L(Lya)= (3.9 ±0.2) x 10 erg s , well i nto the LAB 



categ ory. The Lya emission yields a star- formation rate (SFR) of ~ 36 M0/yr (jOuchi et al. 
2OO9I ): However, this is likely to be an underestimate by a factor of at least a few, given 
that the Lya line is both resonantly scattered by dust, and, for these redshifts, attenuated 
due to absorption by the damping wing of neutral hydrogen in the intergalactic medium. 
The SFR of Himik o is thus likely to be significantly larger than the estimate of 36 Mq/jt. 
Ouchi et al.l (120091 ) used fits of stellar synthesis models to the spectral energy distribution 
over the optical to near-infrared (near-IR) wavelength range to obtain an extremely high 
stellar mass, (0.9 - 5) X 10^° Mq, indicating that Himiko is a massive galaxy. The large 
transverse size of the Lya-emitting region is also consistent with Himiko being an extremely 
massive object, in the context of hierarchical models of structure formation. It is thus an 
excellent candidate for a search for the molecular gas that fuels its high rate of star formation. 



3. Observations and Data Analysis 

Observations of the CO J=2-l line emission in Himiko were made with the 110m 
Green Bank Telescope (GBT) in project AGBTlOB-010, using the Ka-band receiver which 
covers the frequency range 26 — 40 GHz. At the redshift of Himiko, the CO J=2-l line 



^When required, th is paper uses a ACDM cosmology, with Hq = 71 km s ^ Mpc ^, Ha = 0.73, = 0.27 
dSpergel fc et aLll20n7h . 
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[i^rest = 230.538 GHz) is redshifted to 30.354 GHz. The GET beam at this frequency 
has a full- width-at-half- maximum (FWHM) of ~ 25", corresponding to a spatial extent of 
~ 140 kpc at z = 6.6. The GET Autocorrelation Spectrometer was set up to sample two 
800 MHz IF bands at this frequency, overlapping by 200 MHz, with a single polarization 
and 2048 channels per band. The total velocity coverage was ~ 14000 km s~^. Eandpass 
calibration was carried out using sub-reflector nodding with a cycle time of 6 seconds. While 
sub-reflector nodding has high observing overheads, it has been shown to yield excellent 
spectral baselines in earlier GET Ka-band observations, far better than have been obtained 
with standard position-switching. Observations of 3C147 were used to correct the large- 
scale baseline structure in the target spectra and to calibrate the flux density scale, which 
we estimate to be accurate to 10%. Typical system temperatures were 35 K throughout the 
observations. The total on-source time was 4 hours. 

The data were analyzed using the GBTIDii data analysis package, following standard 
procedures. The data were flrst visually inspected and any data affected by correlator fail- 
ures or radio frequency interference were edited out. Next, the narrow expected line width 
(250 km/s, or 25.5 MHz if the CO linewidth is similar to that of the Lya line) relative to the 
size of each 800 MHz IF band, meant that we were able to flt and subtract out a 5th-order 
polynomial baseline from each individual spectrum, without affecting structure on the scales 
of the expected line emission. After flux density and bandpass calibration and baseline sub- 
traction, the data were inspected for any residual baseline structure that might introduce 
spurious signals into the flnal spectrum. Around 10% of the data were found to be affected 
by such artefacts and were edited out. The individual spectral records were than averaged 
together, weighting by the system temperature, to produce the flnal CO spectrum. This was 
then Hanning-smoothed and re-sampled, before further boxcar smoothing and re-sampling 
to resolutions of 50, 100, 150 and 250 km/s. 



4. Results and Discussion 

The two panels of Figure 1 show the flnal CO J=2-l spectra of Himiko at velocity 
resolutions of 100 km s~^ (left panel) and 250 km s^^ (right panel). The root-mean-square 
(RMS) noise is 99 /iJy per 100 km s~^ channel and 68 /iJy per 250 km/s channel. No 
statistically-signiflcant spectral features are visible in either spectrum. 

Our non-detection of CO J=2-l line emission in Himiko places strong constraints on 
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Fig. 1. — Final GBT Ka-band spectra from the z ~ 6.595 LAB, Himiko, covering the 
expected redshifted frequency of the CO J=2-l hne emission; the velocity scale is relative to 
z = 6.595. The spectra in the left and right panels have been smoothed to, and re-sampled 
at, velocity resolutions of 100 km s~^ and 250 km s~^, respectively. The RMS noise on the 
spectra are 99/iJy per 100 km s~^ channel and 68/xJy per 250 km s^^ channel, respectively. 
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the CO line luminosity. The CO line luminosity, L'qq, can be written as 



L'co = 3.25 X 10^ Sco z/, 



2 r)2 

obs 



where fobs is the observii iK frequency (in GHz ), Dl is the luminosity distance (in Mpc), and 
L'qq is in K km s~^ pc^ (ISolomon et al.lll997l ). In the case of a non-detection, 5*00 x = 
3a Av X (in Jy km s~^) gives the 3a upper limit on the integrated flux density in the CO 
J=l-0 line, where a^v is the RMS noise at the velocity resolution AV. We assume a line 
width of AV = 250 km s~^, from the observed Lya line width and simila r to the median 
observed line width in quasar host galaxies at intermediate redshifts (e.g. ICarilli fc Wang 



20061 ). This resuhs in a 3a limit of L'^q < 1-8 x 10^° x (Al^/250)i/2 ^ km s"^ pc^ on the 
CO J=2-l line luminosity. 

Assuming that the CO J=2-l line emission is thermalized, we can use the CO line 
luminosity limit to obtain an upper limit to the cold molecular gas mass in Himiko, as- 
suming an appropriate conversi on factor, Xco, from line luminosity to molecular gas mass 



( ISolomon fc Vanden Boutll2005[ ). This factor can vary between 0.8 and 4.6 Mq (K km s~^ pc^) 
depending on whether the galaxy is an ultraluminous infrared galaxy (ULIRG) undergoing 
starburst activity {Xqo ~ 0.8 Mq (K km s~^ pc~^)~"^; 



Downes fc Solomon 



quies cent object like the Milky Way {Xco ~ 4.6 (K km s~^ pc~^)^"; e.g 
1991 ) . This factor may be even higher in regions of low metallicity gas (e.g. 
2010l : iLeroy et al.ll201ll ). such as those expected to be found within the interstellar medium 
of LAEs in the early Universe. However, massive galaxies like Himiko are unlikely to have 
such low metallicities, especially given its high inferred stellar mass (> 10^° M©). 



1. 



19981 



or a more 



Solomon fc Barrett 



Hughes et al. 



Previous studies of CO J=l-0 line emission in two LAEs at z ~ 6.5 by IWagg et al. 



( I2OO9I ) adopted a low Xco factor, similar to ULIRG v alues, based on the elevated star 
formation rates in these objects. We note that ICombed ( I2OIOI ) argues that the Xco factor 
in these systems should be closer to the Milky Way value. For a ULIRG conversion factor, 
Xco = 0.8 Mq (K km s"^ vc- ^)-\ we obtain the limit M(H2) < 1.4 x 10^° M© {3-a) on the 



molecular gas mass of Himiko. lOuchi et al.l ( 20091 ) obtain a stellar mass of ~ 3.5^2 e ^ 10^° M© 
from a fit to the spectral energy distribution; the molecular gas mass is thus less than the 
stellar mass by a factor of > 2.5. Conversely, assuming a Milky Way conversion factor, 
Xco = 4.6 Mq (K km s^^ pc^^^-i^ yigi^jg ^he limit M(H2) < 8.3 x 10^° M© (3a) on the gas 
mass. 

Our limit on the CO line luminosity in Himiko can be converted to a limit on the total 
far-infrared (FIR) luminosity due to star formation . The best-fit parametrization of the 
relationship between Lfir and L'nn, valid for nearby s piral galaxies and starbursts, is Lpm = 
(1.26 ± 0.05) x logL'co — 0.81 flGao fc SolomonI 12004 ). Assuming that this relationship also 
holds for the interstellar medium of high- 2; LABs, our upper limit to the CO line luminosity 
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yields the upper limit LpiR < 1.3 x 10^^ L© on the total FIR luminosity. Using the relation of 
KennicuttI (Il998[ ) to convert this limit on the FIR luminosity to a limit on the obscured SFR, 
we obtain SFR < 225 Mq yr^^. This is consistent with the lower limit of SFR > 36 Mq yr~^ 
obtained from the Lyman-a line luminosity, even if obscuration/absorption effects cause the 
SFR derived from the Lyman-a line to be underestimated by a factor of a few. 

Searches for CO emission have been carried out i n only two LABs , LABdOS at z ~ 2.656 
and SSA22-LAB01 at ^ ~ 3.09, prior to this work ( lYang et al.l 120121 ). These authors used 
the IRAM 30 m telescope to observe the CO 5-4 and CO 3-2 transitions in LABd05 and 
the Plateau de Bure interferometer to observe the CO 4-3 and CO 3-2 lines in SSA22- 
LABOl. The non- detect ions yielded limits of L'co(5_4) < 1-4 x 10^°(Al^/400)^/2 ^ km s"^ pc^ 
(LABd05) and ^^'^0(3^2) < 1-5 x lO^^AV/my/^ K km s'^ pc^ (SSA22-LAB01), similar to 
our limit on the CO J=2-l line luminosity in the z ~ 6.595 LAB. We note, in passing, 
that similar limits on the CO line luminosity have been obtain ed from GBT CO J=l-0 
studies of two LAEs at z ~ 6.56 and z ~ 6.96 ( IWagg et al.ll2009l ). Thus, at present, there 
is no evidence for si gnificant amounts of molecular gas in h igh- 2: LABs, despite their large 



stellar masses (e.g. lOuchi et al.l l2009l : IColbert et al.ll201ll ). The CO luminosity in LABs 



appears to be significantly lower than that in the l umin ous quasar host galaxy population 
at 2; > 6 (e.g. IWalteret al.ll2003l : IWang et al.ll2010l . l201l[ ). in sub- mm starburst and quasar 
host galaxies at z ^ 2, as wel l as in quiescent B2:K-selected disk gala xies at 2; ~ 1.5 (e.g. 



Solomon fc Vanden BoutI 120051 : IPaddi et al.ll2008l : lAravena et al.ll2010r ). 



Finally, it should be emphasized that the conversion factor from CO line luminosity to 
molecular gas mass is highly uncertain in high- 2; galaxies. A low conversion factor [Xco = 
0.8 Mq (K km s~^ pc^)~^], apphcable to the local ULIRG populatio n, is usually assume d 



for FIR- luminous objects like high- 2; sub-mm galaxies and quasars (e.g. iTacconi et al.ll2006l ). 
Conversely, values closer to ones found in the Milky Way have been obtai ned in star-formin g 



BzK-selected galaxies at 2; ~ 1.5 [Xco 



3.6 ±0.8 Mq (K km s^^ pc^) 



-2. 



Daddi et al. 2010 



While studies of the hi gh- 2: LAB and LAE population hav e typically assumed the ULIRG 



conversion factor (e.g. IWagg et al.l l2009l : lYang et al.ll2012l ). it is quite possible that a high 



conversion factor is applicable to LAEs and LABs at z > 6.5. This may be especially critical 
for the z > 6 LAE population, which appears to consist of relativ ely small galaxies with 



low SFRs, < 10 M(7i yr ^, and low st ellar masses, < 10^ Mq (e.g. iFinkelstein et al.ll2009 



Pentericci et al.ll2009l : lOno et al.ll2010l ). A low CO-to-H2 conversion factor value would imply 
that detections of low- J CO line emission would be extremely challenging in these systems, 
even with new facilities like the Expanded Very Large Array (EVLA) or the Atacama Large 
Millimeter Array (ALMA). For example, the CO J=3-2 line from a 2; = 6.5 LAE with 
SFR=10 Mq yr^^ would be redshifted to a frequency of ~ 46 GHz with a peak flux density 
of Sco ~ 27 /iJy, assuming a line FWHM of AV = 300 km s~^), if we assume Xco = 
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4.6 M0 (K km pc^)~^. Detecting such a line would require hundreds of hours with either 
the EVLA Q-band receivers or the ALMA Band-1 receivers. This i ssue would be exacerbated 
for the z > 6 LAEs if, as it is likely (e.g. iFinkelstein et al.ll201ll ). they are low-met allicity 



factor could be as high as 30 (K km s 



Lerov et al. 


2011 


), 


dels of 


Vallini et al. 



(I2OI2I ) also suggest that the detection of high- J CO line emission in z 
be possible in reasonable ALMA integration times. 



6.6 LAEs will not 



5. Summary 



We present results from a deep GBT search for redshifted CO J=2-l line emission in 
Himiko, a Lyman-a blob at 2; ~ 6.595. We do not detect CO J=2-l line emission, with an 
upper limit of L'(.q < 1.8 x 10^° x (AV^/250)^/^ K km s"^ pc^. T his limit is s imilar to the 



limits on the CO line luminosity obta ined in two LABs a t 2; ~ 3 by lYang et al.l (120121 ) and in 
two Lyman-a emitters at 2; > 6.5 by IWagg et al.l (120091 ). Despite their high stellar masses, 
high- 2; LABs appear to have significantly lower CO line luminosities than luminous quasars 
or sub-mm galaxies at similar redshifts. Our constraint on the CO line luminosity implies 
an upper limit of M(H2) < 1.4 x 10^° Mq on the molecular gas mass, assuming a ULIRG 
conversion factor, Xco = 0.8 Mq (K km s~^ pc^)~^. However, it should be noted that this 
conversion factor is unknown for LABs. A high value of Xco, similar to values in the Milky 
Way or 2; ~ 1.5 B2;K galaxies, would make it extremely difficult to detect CO line emission 
from LABs and LAEs at 2; > 6, even for upcoming telescopes like ALMA and the EVLA. 
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